Fluctuation properties of an effective nonlinear system subject to Poisson noise 
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We study the work fluctuations of a particle, confined to a moving harmonic potential, under 
the influence of friction and external Poissonian shot noise. The asymmetry of the noise induces 
an effective nonlinearity in the potential, which in turn leads to singular features in the work 
distribution. On the basis of an analytic solution we find that the conventional fluctuation theorem 
is violated in this model, even though the distribution exhibits a large deviation form. Furthermore, 
we demonstrate that the interplay of the various time scales leads to critical behaviors, such as a 
negative fluctuation function and a divergence in the work distribution at the singularity. In a certain 
parameter regime large negative work fluctuations are more likely to occur than the corresponding 
positive ones, though the average work is always positive, in agreement with the second law. 

PACS numbers: 05.40.-a, 05.70.Ln, 02.50.-r 
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Fluctuation theorems play a prominent role in the 
study of nonequilibrium fluctuations and have been de- 
rived for a varietyof different systems and thermostat- 
ting mechanisms [H 0, S, 0] ■ In a nonequilibrium steady 
state, which represents the next simplest generalization 
of equilibrium, the so-called conventional fluctuation the- 
orem states that the probability distribution IL(a T ) of 
finding a particular value of a thermodynamic quantity 
a T (usually work or heat) over time r satisfies a certain 
symmetry relation, which, for stochastic systems, can be 
formulated as [H, [j] 



n(-o T ) 



(1) 



where = indicates the behavior for large r. Eq. §T§ 
is a refinement of the second law, in that it quantifies 
the probability of observing temporary second law viola- 
tions (negative a T ) in the steady state. The validity of 
Eq. {!]) is commonly related to the existence of a large 
deviation form of the probability distribution US]. In 
the present letter we further investigate the relationship 
between large deviations and the fluctuation theorem. 
We show, by explicitly calculating the work distribution 
in an effectively non-linear Markovian Langevin model, 
that the conventional fluctuation can be violated even 
though the distribution admits a large deviation form. 
Our model exhibits a range of other striking properties 
such as singularities and critical behavior. These fea- 
tures are all explained on the basis of the time scales in 
the model, as we discuss in detail in the following 

We consider an overdamped particle under the influ- 
ence of a time-dependent force as well as friction and 
external noise. The basic equation of motion for the po- 
sition x(t) of the particle in the laboratory frame reads 



ax(t) = -n(x(t) - vt) + £(t). 



(2) 



Here, the force —K,(x(t) — vt) stems from a parabolic po- 
tential U(x,t) — K,(x — vt) 2 /2 which moves with constant 
velocity v, where either v > or v < 0. The parameter 



k denotes the strength of the potential, a the friction co- 
efficient, and £(t) stochastic noise from the environment. 
Similar models, where £(t) is given by thermal Gaussian 
and external Levy noise respectively have been investi- 
gated in @, 0] ■ If we move the potential for a time pe- 
riod r, a certain amount of work is done on the particle, 
namely 



W T = ~KV / (x(t) - Vt)dt. 

Jo 



(3) 



In the following we are interested in the properties of the 
steady state work fluctuations in the model Eq. ([2]) , when 
the noise is given by a sequence of stochastic kicks 
of variable frequency A and amplitude T. More precisely, 
we consider Poissonian shot noise of the form [81 



:(t) = J2 T kS(t-t k ), 



(4) 



fc=i 



where n t , the number of kicks in time t, is determined 
by the Poisson counting process P(n,t) = (\t) n e~ xt /n\. 
The parameter A denotes the average number of kicks 
per unit time (rate of kicks) so that there are Xt kicks 
occurring in the time interval [0, t]. When a kick occurs, 
the amplitude T k is sampled randomly from a distribu- 
tion p(r). In the following we focus on one-sided shot 
noise and assume an exponential distribution of ampli- 
tudes p(T) — Fq e~ r / r °. Noise specified according to 
Eq. (|2]) is characterized by delta-correlated cumulants 
and mean (z(t)) = XTo [9(. It is convenient to use a 
noise with a zero mean so that we consider in Eq. ([2]) 



m = z(t) - Ar . 



(5) 



Using this form of the noise and changing to coordi- 
nates in a comoving frame, y(t) = x(t) — vt, the equation 
of motion © reads 



m = — vit) 



AFn \ 1 , . , * 

— ~ + -*(t), (6) 
a I a 



2 



and the work Eq. ([3]) is given as 



the Poissonian shot noise: A 



■i<a I y(t)dt. 
lo 



(7) 



In Eq. ([5]) r r = a I n denotes the characteristic relaxation 
time of the oscillator. It is crucial that Poissonian shot 
noise gives rise to two additional time scales in our model, 
namely the characteristic time scale of the fluctuations 
t\ = A -1 , which is the mean waiting time between two 
successive kicks, and the time scale 



To 



(8) 



which relates the mean amplitude of fluctuations and the 
friction due to the driving. We will see below that the 
three time scales r r , t\, and r p determine the properties 
of the work fluctuations in our model, and that critical 
behavior appears due to the interplay of these times. 

From the Langevin equation ([6|) we can infer two im- 
portant properties of the model. Firstly, upon averaging 
of Eq. ([6]) we obtain the mean position in the steady 
state (y(t)) = —vr r and from Eq. |(7J) the mean work 
(W T ) — av 2 r, which is always positive in agreement with 
the second law. Secondly, we find that there exists a 
minimal value y* of the position coordinate. This can 
be seen if we solve Eq. © without the stochastic term 
z(t), which yields y* = — (v + XTo/a)T r in the long time 
limit. The important observation is that the influence of 
the noise only provides kicks in the positive direction (cf. 
Eq. (@|) so that y* is a cut-off, i.e., no positions below 
y* can be reached. Since the work rate is proportional to 
the position, we likewise obtain from Eq. ([7|) an extremal 
value of the work, namely 



W* T = (W T ) (1 ± Xt p ), 



(9) 



where the + sign corresponds to v > and the — sign to 
v < 0. It is important to note that the spatial asymme- 
try of the noise induces a qualitative different behavior 
of the work fluctuations depending on the sign of v. In 
the case v > the value of W* is always positive and 
corresponds to the maximum work done on the system 
in time r. For v < on the other hand, W* is the mini- 
mum work value. In that case W* can be either positive 
or negative (cf. Eq. ([7])). Its sign is determined by the 
two time scales t\ and t p . From Eq. (J9j) we see that W* 
is positive if t\ > t p and no negative work fluctuations 
can occur. Furthermore, W* — if t p — t\ and W* < 
if t\ < T p . The position cut-off y* represents an infinite 
barrier in the potential so that the noise induces an ef- 
fective non-linearity in the potential. A work cut-off can 
also be observed in a Brownian particle model, where the 
moving potential is give n as a non-linear potential of the 
Lennard- Jones- type [ill ]. Furthermore, due to these cut- 
offs, the distributions of position and work are generally 
non-Gaussian unless one considers the Gaussian limit of 



AFn = const. In this limit y 



oo and To — > with 
o — uuimi. ±11 umo mi™ y —r — oo and W* — > ±oo. 
In order to determine the work distribution for arbi- 
trary A and T values we first derive an exact expression 
for the characteristic function of the work using a theo- 
rem on generalized Ornstein-Uhlenbeck processes 0, E3] ■ 
This theorem states that the characteristic functional of 
the process y(t) of Eq. ©, defined as 

G v [h(t)] = (exp{iJ h(t)y(t)dtj\ , (10) 
is determined from the noise functional via 



G y [h{t)]=e i y° k °G i ,[k{t)], 



(11) 



where G^ denotes the characteristic functional of the 
noise plus the drift terms in Eq. ([6]): £'(t) = £(t)/a — v. 
Furthermore, the functions fco and k(t) are related to the 
test function h(t) according to k(t) = h(s)e^~ s ^ Tr ds 
and fco = J °° h(s)e~ s / Tr ds [T3|. It is then important to 
note that choosing the particular test function h{t) = 
-qnvQ(T - t) in Eq. (JTOj) leads to 

G y [h{t)\ = <^exp j^g^ y{t)Atyj = G Wr (q), (12) 

i.e., the characteristic functional for h reduces to the 
characteristic function of the work Gw T (o) = {e lqWT ) due 
to Eq. © 0. In our case the noise functional G^[k(t)] 
is given by 

G e [k{f)\ - e- i{v+Xr ° /a) ^ mdt G z [k(t)/a], (13) 

where G z denotes the characteristic functional of the 
Poissonian shot noise z(t) whose exact analytical expres- 
sion is known [8fl: 



Gz[k{t)] = exp <^ A 



((• 



,irk(t) 



l)dt 



(14) 



Using this noise functional together with the functions 
fco and k{t) obtained with h(t) yields, after some manip- 
ulation, the characteristic work function 



G Wt oc f 1 + iqT v(l - e- T ' Tr ) 



(15) 



exp <j -iqavy* — + Ar ( 1 1 

T r \ 1 + iql v 



with a pole at q — i/iTaav). In the derivation of Eq. (jTB]) 
we have made use of the exponential distribution of am- 
plitudes r and furthermore chosen an initial condition yo 
sampled from the noncquilibrium steady state distribu- 
tion, which can be found by solving the Fokkcr-Planck 
equation associated with the Langevin equation (JSJ) 9]. 
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For our discussion of the work fluctuations we intro- 
duce the scaled dimensionless work value p, defined by 
p = W r j (W T ). The distribution of p is obtained from 
the inverse Fourier-transform of G\y T 



n T (p) = ^v- 



2tt 



G WT {q)e~ 



iqpav t 



dq, (16) 



where N denotes the normalization constant. To our 
knowledge there is no exact result for the inverse Fourier 
transform Eq. (|16p . However, for large r the integral will 
be dominated by its saddle-point and can then be analyt- 
ically obtained using the method of steepest descent [l2j . 
Neglecting terms of order t^ 1 / 2 then yields the following 
result for the distribution II T (p): 



n T (p) 



4tt \p* - 1| 



x exp < 

ta 



jp* 


-p 


p* 




p* - 


-p 


p* - 


- 1 



- 1 



(17) 



where p* denotes the rescaled extremal value of work: 
p* = W*/ (W T ) = 1 ± Tp j T\ (cf. Eq. ©). Here, as well 
as in the rest of this letter, we have expressed A as 1 
in order to emphasize the crucial role of the time scales. 
Eq. (fl7| shows that the distribution H T (p) is completely 
specified by the times r r , t\, t p , and r. Importantly, 
the square root \/ (p* — p)/ {p* — 1) is always real, since 
both p* — p and p* — 1 are either positive (v > 0) or 
negative (v < Q). For \p* — p\ ^> the tail decays faster 
than exponential but becomes exponential for very large 
r, where Tl T (j>) exhibits the large deviation form 



n T b) = e 



with rate function 



Hp) 



1 



P ~P 
p* - 1 



- 1 



(18) 



(19) 



Furthermore, one notices two different singularities ap- 
pearing in Eq. (JXTJ) . Firstly, the derivative of n T (p) di- 
verges for p — > p* as n'(p) cx \p* — p\ in leading order. 
This means that the approach of n r (p) to the cut-off has 
a vertical slope (see inset Fig. QJ . Secondly, one notices 
that n T (j?) itself diverges for p — > p* if -tv/ia < 3/2. This 
divergence for small t t /t\ at p = p* can be understood by 
considering the nature of the work and position cut-offs 
in our model. If T r /r\ is too small the system relaxes 'too 
quickly' in between the stochastic kicks and thus spends 
most of its time at the position that it would assume 
without noise, i.e., at y*. Consequently, the particle will 
predominantly acquire work W* over time r leading to a 
divergence in the work distribution at p — p* (see inset 
Fig.QJ. The condition t t /t\ < 3/2 can also be expressed 
in terms of a critical friction coefficient a* = |kta, so 
that the divergence of H T (p) appears when a < a* . 




FIG. 1: Comparison of the analytic saddle-point approxima- 
tion Eq. (|17[l with a numerical inverse Fourier transform of 
Gw T and results from a direct simulation of the equation of 
motion (|6]). Parameter values: r = 10r r , r r = 1, v = 1, 
t\ — 0.2, To = 0.5. Inset: The approach to the cut-off given 
by Eq. (|17[) . For t t /t\ > 3/2 the distribution goes to zero 
(black line), while it diverges when t t /t\ < 3/2 (blue line). 
Parameter values: r — 10, r r = 1, t v = 1. 



We find that Eq. (jXTJ) yields an excellent approximation 
of the distribution IL T (p) at least for r > 10r,.. This 
becomes evident in Fig. [TJ where we compare Eq. (|17|) 
with a numerical inverse Fourier transform of Gw T an d 
also with results from a direct simulation of the equation 
of motion ([6]) using a Poissonian increment method [l3T | . 

In the asymptotic regime r — > oo the work distribu- 
tion Eq. (|17p is dominated by the large deviation form 
Eq. (fT8"|) . In order to further discuss the fluctuation prop- 
erties of work we consider the fluctuation function 



fr( P ) 



1 



n r (p) 

(W T ) "* Ilr(-P)' 



In 



(20) 



The conventional fluctuation theorem then predicts that 
linv^oo f T (p) = p. From Eq. (fT5|) we obtain instead 



lim f T (p) 

T — >00 

= 2p | 2(p*-l) 
vT Q vT q 



p -p 



p" 



1 



p + p 



p" 



(21) 



defined on the interval [— p*,p*\- This restriction on the 
range of p- values is similar to deterministic Anosov sys- 
tems, where the phase-space is bounded 0]. We see that 
the conventional fluctuation theorem is violated in our 
model, even though we have identified a large deviation 
form of the distribution. We characterize the behavior of 
the work fluctuations for v > and v < separately. 

(i) For u>0we find that Eq. (f2"Tj) has the zeros p\ = 
and P2 = 2-^/p* — 1 (due to the symmetry of f T (j>) we ne- 
glect the negative root) . The zero p2 becomes significant 
when p* > 2, because then p* > p2 and p can assume 
values in the interval [j»2,P*]- The crucial observation is 
that linv^oo f T (p) from Eq. (pH]) becomes negative for 
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FIG. 2: The fluctuation function g(p). At the cut-off g'(p) 
diverges while g(p) assumes a finite value, a) v > and four 
different p* values. For p* > 2 we observe that g(p) becomes 
negative in the interval p G [pojP*]- b) v < and four different 
p* values. Here, g(p) only has a zero at p = 0. 



p G [p2,P*], if P* > 2. This is evident in Fig. [2^), where 
we plot the dimensionless rescaled fluctuation function 
g(p) = lim T ^oo \v\Tof T (p), which depends only on p*. 
There exists therefore a parameter regime in which neg- 
ative fluctuations of a certain magnitude are more likely 
to occur than corresponding positive ones. In fact, since 
p* = 1 + t p /t\ we find that p* > 2 if r p > T\. Using the 
definition of t p , Eq. ([5]). we can identify a corresponding 
critical velocity v* = T /(aT\) so that for < v < v* the 
fluctuation function is negative. 

This surprising property originates from the strongly 
asymmetric tails of the work distribution H T (p): The 
negative tail decays exponentially while the positive tail 
is bounded by the cut-off at p*. It is important to note 
that despite the existence of a negative regime of the 
fluctuation function f T {p), the second law is never vi- 
olated: the mean value of the work is always positive, 
(W T ) = av 2 r. For p — > p* the derivative g'(p) diverges 
like (p* — p)^ 1 / 2 , while g{p) remains finite (see Fig. [2k))- 

(ii) For v < we have to distinguish two regimes. 
Firstly, for t\ > t p no negative work fluctuations can 
occur (cf. Eq. ([9])) and therefore the fluctuation function 
f T {p) can not be defined in this parameter regime. Sec- 



ondly, for t\ < T p we have p* < and we can discuss 
g(p) in the interval [0, —p*]. In this case g(p) is only zero 
at p = 0. As in case (i), the derivative g'(p) diverges for 
p — ► p* while g(p) remains finite (see Fig. HJa). 

We note that our theory could be adapted to an ex- 
periment similar to that of Mahadevan et al where 



a lubricated rod of a hydrogel sliding on a soft vibrat- 
ing substrate is considered as a model for biomimetic 
ratcheting motion. Instead of the purely oscillatory vi- 
brations of 14J one could induce asymmetric Poissonian 
shot noise, which could lead to work fluctuations with 
similar features as those presented in this letter. 

In summary, we have investigated the work fluctua- 
tions of a particle, confined to a moving harmonic po- 
tential, under the influence of friction and Poissonian 
shot noise. The one-sidedness of the fluctuations leads 
to a cut-off in the work distribution, which is therefore 
strongly non-Gaussian. We have shown that the work 
distribution exhibits a large deviation form but that, 
nevertheless, the conventional fluctuation theorem is vi- 
olated. 
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